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Abstract 

A recently found [3] 2-index, symmetric, trace-free, divergence-free tensor is 
£^ ■ introduced for arbitrary source-free electromagnetic fields. The tensor can be con- 

structed for any test Maxwell field in Einstein spaces (including proper vacuum), 
*q ' and more importantly for any Einstein-Maxwell spacetime. The tensor is explicitly 

given and analyzed in some special situations, such as general null electromagnetic 
5h . fields, Reissner-Nordstrom solution, or classical electrodynamics. We present an 

explicit example where the conserved currents derived from the energy-momentum 



X 



tensor using symmetries are trivial, but those derived from the new tensor are not. 

1 Introduction: electromagnetic fields 

Let (V,g) be a 4-dimensional spacetime with metric g (signature -,+,+,+). An electro- 
magnetic field is any 2-form F satisfying the Maxwell equations (e.g. fTf\ !2Hj) 

dF = 0, 5F = -j ^ V [fl F up] = 0, V a F° v =f (1) 

where j is the electromagnetic current 4-vector. As usual, one can define the dual of F 
and the associated complex self-dual 2-form J- by means of 

F^\^ P rF p \ T^ l -(F-iF) 

where rj is the canonical volument element 4- form. A duality rotation (e.g. is any 

transformation of the type 

T 1 = e lQ T. 



*Dedicated to Professor Jose A. Azcarraga on the occasion of his 60th birthday. 
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The source-free Maxwell equations — that is, with j = — , which can be written simply 

as V ' a T av — 0, remain clearly invariant under duality rotations. 

The energy-momentum tensor of F can be written in any of the following forms 

I / * * \ 

F^pF v ^ — (jy^-F p a F^ — yF^pFp^ -\- FppFu ^ J 2,J~ apF v p (2) 

where an overbar denotes complex conjugation. Obviously, is symmetric, traceless, 
and invariant against duality rotations. Furthermore, it satisfies the dominant energy 
condition, that is to say 

Too > \Tfj, v \ in all orthonormal bases 
or equivalently, see e.g. [13 EI] 

T(u, v)>0 

for all future pointing (ergo non-spacelike) vectors u and v. Hence, if the energy density 
with respect to an observer vanishes, then the whole energy-momentum tensor is zero, and 
the electromagnetic field is absent. Another important relation satisfied by this tensor is 
the algebraic Rainich identity [26 | l24 l l25| l82 l l6j 

Tp P T% = \{T ap T^)g, v (3) 

which fully characterizes the explicit form (J2J) if T M M = is taken into account. 
A trivial calculation using (0) and (J2J) leads to 

V M T^ = F up j p whence: j = V ^ = . (4) 

This encodes the "covariant" conservation of energy and momentum for electromagnetic 
fields in the absence of charges and electric currents. It also implies that, for any conformal 
Killing vector!; the corresponding Killing current J is divergence-free: 

J»(£) = T^tu => Vp3» = (5) 

providing conserved quantities via Gauss' theorem (THl I2H]- see e.g. [TBI EI] for some 
explicit examples. On the other hand, if there are electric charges and currents on the 
spacetime, one can still restore conservation of the total energy and momentum by adding 
the corresponding energy-momentum tensors of those sources. A typical example can be 
found in ^Zj for classical electrodynamics. 



2 The main result: H^ v 

For any electromagnetic field, there exists another tensor quadratic in F and with the 
same properties as Tp U except for the dominant energy condition. This tensor can be 
written in any of the following equivalent forms 

= V T Fp P V T F/ - \gp U V T Fp a V T F^ = 

= \ (V T F W V T F/ + V T F, P V r F v p ) = 2 V T ^ P V T T v p (6) 
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from where immediately follows that H^ v is symmetric, traceless, and invariant under 
duality rotations (with constant 9 !). More importantly, it is also divergence-free if the 
source-free Maxwell equations hold 

V^iP" = V T F vp V T j p whence: j = V^tf"" = . (7) 

This is valid (i) if the electromagnetic field is "test" , that is to say, it does not contribute 
to the righthand side of the Einstein field equations, so that the spacetime is vacuum 
with a possible cosmological constant. And (ii) in the full non-linear theory such that the 
Einstein-Maxwell field equations with a possible cosmological constant are satisfied. The 
proof of these results is far from obvious using tensor calculations. Actually, its derivation 
was achieved using spinors in [S] . For a full proof involving only tensor techniques one can 
consult [12] (see also [TT] , where related results for a "Lanczos potential" L^ UT in place of 
Vt-F^ were found). 

Analogously to the case of T^ u , using (JZJ one can construct divergence- free vector fields 
associated to any conformal Killing vector £ by means of 

r(0 = = o (8) 

providing new conserved quantities in any Einstein-Maxwell spacetime (or any Einstein 
space with a test F) having such £. 

Note finally that, generically, the dominant property (i.e. that H(u,v) > for all 
future-pointing u, v) need not be satisfied. 

An obvious question arises: why was unknown until very recently? One reason 
may be, as previously stated, that the proof of (JZJ) involves difficult calculations. Another 
one may be the circumstance that its physical units are not those of energy density. 
Actually, has units of "superenergy" in analogy with the well-known Bel-Robinson 
tensor, see e.g. j2 El HI El E3 EI], which means units of energy density per unit surface. 
One can further prove (see [3]) that in the unique trace of the 4- index superenergy 
tensor of the electromagnetic field, as given originally by Chevreton jH] or in the unified 
general treatment of [HU EH]- Thus, if the new tensor has any physical relevance at all, 
this could help in understanding the concept of gravitational superenergy as well. 

3 Uniqueness of H^ u 

The tensor © is unique and solely constructable for the source-free electromagnetic field. 
I would like to stress these properties by presenting a list of more detailed points: 

1 In Special Relativity, that is to say, in flat spacetime, there exist many tensors which 
are quadratic in F and with the previous properties. One example is the so-called zilch 
tensor [2T1 [THJ , but there are many other. For a full list, see [HUH]. Several examples 
pertinent here are 

H$ = V T1 . . . V Tn F^V- . . . V-F/ - \g, v V Tl . . . V Tn F p(T V^ . . . V*F>° 
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for any natural number n. Observe that H P J = T pv and H P J = H^. Obviously, these 
tensors are traceless and symmetric, and furthermore they are divergence-free for any 
n if the electromagnetic field has no sources. This is obvious from the commutation of 
the covariant derivatives in flat spacetime. However, the only tensors of the previous list 
which keep the divergence-free property in curved spacetimes are T pu and H^. 

2 The divergence-free property of H pv in generic spacetimes holds in four dimensions 
exclusively. Thus, (J7J) is not valid, not even for test electromagnetic fields, in higher 
dimensions. 

3 The construction of H pi/ is not translatable to other physical fields. Take, for instance, 
a minimally coupled massless scalar field <fi. The tensor analogous to H pu can be seen to 
be 

= V P V^ V p V u <P - X -g^ V p V CT VV>. (9) 

This can be derived in several independent ways, such as by taking the trace of the 4- 
index superenergy tensor of the scalar field jSH EDI EH], or in direct analogy with the 
expressions © and ^ if one recalls the energy-momentum tensor of 0: T pv = V M 0V^0 — 
I^Vo-^V "^. Nevertheless, in non-flat spacetimes, and despite assuming the massless 
Klein-Gordon equation D0 = 0, D pu is not divergence-free. This happens even for test 
scalar fields, that is to say, in Ricci-flat spacetimes. 

4 The existence of H pv in the full Einstein- Maxwell theory, specially as given in its 
form (jHJ) not involving any curvature terms — just first derivatives of F- may seem a little 
surprising. For source-free test electromagnetic fields in Einstein spaces — including the 
Ricci-flat case — , the existence of can be understood as a "second derivative" of the 
energy-momentum tensor. To see this, let us start with the simplest case of flat spacetimes, 
where clearly 2H pv = DT^ — if the source-free Maxwell equations ((TJ) hold — due to the 
commutativity of the covariant derivatives. In curved spacetimes things are not so simple. 
However, in Einstein spaces, there is a D'Alembertian-like operator introduced by 
Lichnerowicz [20] which commutes with the divergence. Thus, after some manipulations 
and using some tensor identities one can actually prove that 2H pu = D L T pv in Einstein 
spaces. This result has been pointed out in for the Ricci-flat case. 

One can thus ask whether a similar but more involved calculation permits a relation 
between the energy-momentum tensor and the new tensor via a second order operator 
in the full non-linear Einstein- Maxwell case (allowing for an additional cosmological con- 
stant). This has been solved by Edgar [TO] recently in an elegant manner. He shows that 
one can define in general an operator acting on divergence-free rank-2 tensor fields such 
that the result is divergence-free too. Nonetheless, this operator involves curvature terms 
in addition to second order derivatives. Thus, the path to prove the relation of this with 
H pu is still complicated. As a matter of fact, in order to rederive the expression © one 
has to (i) restrict the discussion to 4 dimensions, (ii) use the Einstein field equations, 
(iii) then utilize the source-free Maxwell equations (0) explicitly, (iv) take into account 
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the special Rainich identity 0, and (v) manipulate the result employing a 4-dimensional 
identity fUE] involving the curvature tensor and the electromagnetic field 1 . 

One could argue, as for instance in that this deprives H^ v of relevant "independent 
content" from (J2J). This is not clear to me, though, for there are cases in which the 
divergence-free currents (JSJ) constructed from the energy-momentum tensor are trivial 
while those built from the new tensor, as in (jHJ), are not. An explicit example is presented 
in ^4.1. 21 The point here is that the commutativity of the divergence and the second 
order operator does not immediatey translates to the currents, as they involve (conformal) 
Killing vectors. Thus, the curved space derivatives of the infinitesimal symmetry are also 
involved here. 

4 Applications and explicit examples 

In this final section we present some explicit expressions for the new tensor H^ v in some 
cases of physical interest. The results may shed some light into its interpretation and its 
potential usefulness. 

4.1 Null electromagnetic fields 

Null electromagnetic fields can be characterized by any one of the following equivalent 
conditions: 

T^T^ = or equivalents F^F^ =F^ V F^ = 0, 

3/: FF^ = F F^= (then necessarily tH^ = 0), 

F = £Ap, F= £ A q (10) 

where the unit spacelike vectors p and q are orthogonal to each other and orthogonal to 
£, the latter being intrinsically defined (up to a proportionality factor) by the null field 
and sometimes called "wave vector" — it defines the null direction of propagation of the 
electromagnetic field — . It must be noted that p and q are not intrinsically defined by F 
and are subject to redefinitions of the form p — > p + a£, q — > q + b£ for arbitrary functions 
a and b. In other words, only the direction of the wave vector i is intrinsically defined and 
therefore, its orthogonal spacelike 2-plane, which is spanned by p and q and sometimes 
called "polarization plane" , is also intrinsic to the null F. 

Due to the above remarks, the energy momentum tensor should only depend on the 
wave propagation direction, and this is certainly so, as is well known, because using (J2J) 
and (fTUj) one immediately gets 

T = £ £ 

1 Points (i) to (iv) could be analogously followed for the scalar field without problems — substituting 
the Klein-Gordon equation for the Maxwell equations — , because there is an identity similar to @ in 
that case too, see e.g.jH]. However, the final point (v) has no translation to the scalar field, and thus 
an expression without curvature terms, such as I©, is impossible for divergence- free tensors in this case. 
This holds also for test scalar fields, as remarked in point 3. 
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Hence, the full Einstein-Maxwell system of equations is given by (R^ is the Ricci tensor, 
units with 8ttG = c = 1, and the cosmological constant A is kept) 

together with (JTJ) appropriately specialized on using (fTUj) . An old result due to Mariot and 
Robinson [221 EH] states that this last part requires a null vector field £ which is geodesic 
and shear- free (see e.g. for these definitions), i.e. 

rv M r oc r, v^vt = \iy^f ■ 

But then a series of results eventually summarized in a paper by Goldberg and Sachs ( [T3j 
and references therein) imply that actually the spacetime must be "algebraically special" 
in the sense that the Weyl tensor C a p^ v has a repeated principal direction given by £, see 
e.g. f32] P- 88, that is to say 

For the Petrov classification and the nomenclature about principal null directions, see e.g. 
|U EEl E2! • Using the Newman-Penrose notation (e.g. [25, 32J) these can be written as 

= ^i = 0, a = k = 

in a null tetrad with i as first null vector. The first two here express that £ is a multiple 
principal null direction of the Weyl tensor, and the second pair that it is geodesic and 
shear-free. 

Let us turn our attention to for these null fields 2 . If one plugs (fTUjl into ©, a naive 
calculation leads to expressions involving derivatives of p and/or q. But this is undesirable, 
as we know that only £ is intrinsically defined, so that an expression involving only £ must 
exist. This can certainly be achieved, after a very long and not easy calculation where all 
the previous results and other identitites must be used. The result is 

= V p [i {l N p £ v) - £ p V (lx £ v) - £ {t N v) i p ] . (11) 

This is the general expression of H^ v for null electromagnetic fields. Several alternative 
but equivalent formulae can be given, such as 

= ~n(yy + r&c^ - ^a 

which may also be useful in some particular situations. Nonetheless, this last expres- 
sion depends explicitly on the curvature of the spacetime, and thus (jllj) is preferable in 
principle, and physically more interesting. 

Some particular relevant examples of null electromagnetic fields are examined next. 

2 An analysis of the 4- index Chevreton tensor for null electromagnetic fields can be found in |34|. 
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4.1.1 Maxwell pp- waves 

The plane-fronted waves with parallel rays (pp-waves in short) are defined as those space- 
times containing a constant null vector field £, see e.g. |32] anc ^ references therein, that 
is 

V M 4 = o . 

They contain a large subclass of solutions depending on two arbitrary functions which 
satisfy the Einstein-Maxwell equations. As is obvious from the previous formula and (|TT| 
one has for all these solutions H uv = 0. Thus, the new tensor is identically zero for the 
usual plane waves, including here the classical case of plane waves in Special Relativity. 

4.1.2 Robinson- Trautman type-D solution for a null Maxwell field 



In 1962, Robinson and Trautman [29] found, among many other algebraically special 
solutions, a particular spacetime which solves the Einstein-Maxwell system with a pos- 
sible A for a null electromagnetic field. The solution is of Petrov type D and has a 
3-parameter group of isometries acting transitively on spacelike 2-surfaces. In local coor- 
dinates {u,r, x, y}, the line-element takes the explicit form 

ds 2 = r 2 {dx 2 + dy 2 ) - 2dudr + + ^r 2 ) du 2 

where m(u) — the "mass" — is an arbitrary but non- increasing function of the retarded 
time u: m < (this solution is reminiscent of the more popular Vaidya's spherically 
symmetric radiating solution). The null electromagnetic field is given (up to duality 
rotations) by 

F = V— 2m du A dx 



and its wave vector field is 



du 



r 

so that the only non-zero component of the energy-momemtum tensor is 

2m 

uu ^2 ' 

Using (|TT|) . a direct computation leads in this spacetime to 



mr r A 6 



where H ll v takes the simple form 

\H» V ) =-— xdiag(-l,-l,l,i; 
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which is a structure analogous to that of typical "Coulombian" fields (see subsection 
below). Thus, the form of splits as the sum of a "Coulombian-like" term, due to 
a "charge" represented by ^—2m (the "mass loss" or radiation rate) plus a typical null 
term with contributions from m, m/rh and the cosmo logical constant. 
Using now the three Killing vectors 

& = {d x ,d y ,yd x - xd y } 

one can construct conserved currents in the way given in (JHJ) and (JSJ). It is remarkable 
that, as is obvious, all the currents (0) constructed from the energy-momentum tensor 
are identically vanishing, 

J»{&) = 0, 

whence the corresponding conserved quantities are trivial. On the other hand, the 
divergence-free currents (jHJ) constructed from the new tensor are non- vanishing in general. 
For any one of the three previous & one has 

AS) = ~l , 

leading to non-trivial conserved quantities, as can be easily checked. Thus, in this par- 
ticular spacetime, the new tensor provides useful information, and conservation laws, 
involving the physically relevant magnitude fa. 

4.2 Non-null electromagnetic fields 

In this case, it is not possible to give a general formula such as (jllj) because there are two 
null directions defined by F, and they need not be shear-free nor multiple null directions 
of the Weyl tensor. We consider briefly two simple examples. 

4.2.1 The Reissner-Nordstrom solution 

The more general spherically symmetric solution of the Einstein-Maxwell equations is 
given by the following line-element in typical spherical coordinates {t, r, 8, ip} 

ds 2 = -(l- — + < Q\dt 2 +(l- — + < Q\ 1 dr 2 + r 2 (d9 2 + sin 2 9 dp 2 ) 

where m and q are arbitrary constants representing the total mass and electric charge 
of the body creating the gravitational field. Only the exterior asymptotically flat region 
defined by 1 — 2m/ r + q 2 /r 2 > will be considered. This restricts the range of r to r > 
if m 2 < q 2 , or to r > r + = m + ^m? — q 2 if m 2 > q 2 . Then, the metric is static (d t is a 
hypersurface-orthogonal timelike Killing vector) and spherically symmetric. 
The electromagnetic field is given by 

F = -L dt A dr 
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so that its energy-momentum tensor in the coordinate basis is 

(T^) = ^xdiag(-1, -1,1,1). 
A straightforward calculation leads to 

(JP„) = ^ ( 1 - — + 4) X dia 8(- 3 ' 2 > 2 ) • 

Notice the physical units. This particular tensor satisfies the dominant property as well. 
The conserved currents that can be constructed from here are similar to that of the 
energy-momentum tensor, see also [18j. 



4.2.2 The "Bertotti-Robinson" solution 

As early as 1917, Levi-Civita JH] gave the only Einstein-Maxwell spacetime which is 
homogeneous with homogeneous electromagnetic field. However, this line-element is com- 
monly known as the "Bertotti-Robinson" solution, as it was rediscovered in [27] and, with 
the cosmological constant, in [7j. Locally, it can be given in adequate coordinates by 

ds 2 = —dt 2 + dx 2 + cos 2 (ax)dy 2 + cos 2 (at)dz 2 

and the electromagnetic field is 

F = V2a cos(at)dt Adz ==► VF = 

so that H^y is identically zero. The energy-momentum tensor is in this case 

(T^) =2a 2 x diag(-l, -1,1,1). 



4.3 Classical electrodynamics 

In flat spacetime of Special Relativity, and using an inertial reference system, one can 
define the associated electric E and magnetic B fields. Then, using dots for time deriva- 
tives and t i for derivatives with respect to the spacelike rectilinear coordinates x l , one can 
obtain (0 is the time-component and • is the usual scalar product of 4-vectors) 



Hoi 



E ■ E + Ef E' 1 + B ■ B + Hi- H> 



■□(£ • E + B ■ B) 



E x B + IjX S>< = B(E x B) 
2 3 
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